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ABSTRAGT 
Elastic buckling loads are obtained for a family of tapered 
columns using the Vianello iteration method as refined by Newmark. 
The family of columns studied by Dinnik is generalized to include 
columns not symmetric about mid-length. An extensive table of 
buckling factors is compiled with the aid of a high speed digital 


computer. Some of Dinnik's results are found to be in error. 
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l. Introduction 

For an initially-straight slender column which obeys Hooke's 
law, the Buckling load.is defined to be the concentric compressive 
load under which the straight configuration becomes neutrally 
stable. The buckling load is a function of geometry, material, 


and end conditions, and it takes the form: 


sE 
Р = 2 I 1 





where P is the buckling load, m is the buckling factor, E is the 
modulus of elasticity of the column material, I is the minimum 
centroidal moment of inertia of the column cross sectional area, 
and 4 is the column length. 

Columns with varying cross-sections are often used when the 
prismatic column is inadequate. The varying cross-section 
column is used in many applications in order to decrease weight 
or provide additional strength. Early studies of the effect of 
varying cross-sections on buckling loads were done by Dinnik (1) 
and by Bairstow and Stedman (2). Dinnik presented a table for 
columns which tapered symmetrically about a central prismatic 
section. He based his analytical calculations on the assumption 
of various power laws for the stiffness in the end portions of 
the columns. 

Bairstow and Stedman presented a numerical method of solving 
For buckling loads in columns with varying cross-sections. This 
method was later improved by Ratzersdorfer (3). Vianello (4) 
presented a graphical method of solution for buckling loads, and 


Morley (5) introduced numerical integration into Vianello's 


15 


process. Newmark (6) introduced a number of refinements to this 
process and summarized the entire procedure. Approximate 
solutions based on the energy method were given by Bleich (7), 
and Boyd (8) conducted both theoretical and experimental studies 
on tapered columns. 

The present study uses a numerical method of solution to 
Find elastic buckling loads for columns of varying cross-section. 
Specifically, the columns are tapered, hinge-ended, slender, 
initially-straight, and concentrically loaded. А large family 
of columns is studied, and Dinnik's columns, tapered symmetrically 
about mid-length, are members of this family. Dinnik's data are 
checked and augmented by obtaining additional points and by pro- 
viding four significant figure results. The major portion of 
the results, however, provides buckling loads for columns which 
are not symmetric about midlength. These latter results are 
believed to be neu. 

The work can best be described by dividing it into three 
phases: selection of a numerical method of solution; adaptation 
and use of the selected method in solving for the buckling loads; 


and discussion and verification of results. 
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Ze Selection of a Numerical Method 

The selection of a numerical method involved the investiga- 
tion of the Bairstow-Stedman, the Vianello, the Rayleigh (9), 
and the Rayleigh-Ritz (10) methods. Each of these methods is 
designed to find the buckling factor for a varying cross-section 
Eg umm C11). 

In the discussion of the four methods, it is assumed that 


I the moment of inertia at a longitudinal distance x, is known 


x? 
at all points. It is further assumed that the boundary conditions 
are those of a hinged-end column: 
πλου = ডি) 
ក្រយា (272) 
Both the Bairstow-5tedman and the Vianello methods utilize 


the differential equation for the deflection curve of a slender 


column under a compressive load and with small deflections: 


ES zy សយ. 5. [2 5) 


where y is the lateral deflection at the longitudinal distance 
x. It is convenient to modify this equation by a form of 


ВАША 71017012150 





wf 
រ 
Where I, is chosen to be the maximum moment of inertia in the 


2 


column. Equation (2.3) then becomes: 


MS 


22 "р (Ф) у= О 2215 


[55 


In the Bairstow-Stedman method, a value for m is assumed, 
ahd after the quantity т12/1, is calculated, the following 


approximate equations are used in a numerical integration: 


χι - es n A у: ‹ (25) 
% _ — 7, 
Yes 52% 5 У: (2.6) 


Уг... Е ек Б У: STD 


The primes indicate differentiation with respect to x, and y'= 

у! > The column is divided into k increments of equal 

length X, or{=kX. The subscript i indicates the longitudinal 

position and is equal to kxf# , i.e., when i = 0, x = D and when 

i =k, x =. The initial entries Vo and ус are equal to zero in 

accordance with the first boundary condition and equation (2.6). 

In order to begin the integration an initial value for yj at x = s 

is arbitrarily assumed. Then the equations (2.5), (2.6), and (2.7) 

are successively employed until the integration has proceeded 

along the entire length of the column. If, at the completion of 

integration, the deflection corresponds to the second boundary 

GmNEDLLIOn, l.es; ns 0, the correct value of m uas assumed. If 

the end deflection does not meet the boundary condition, new values 

ОГ m must be assumed: until the second boundary condition is met. 
In the Vianello method, m initial deflection curve is 

assumed. Then, using equations (2.6), (2.7), and (2.5), respec- 

tively, the curve yk*/m is obtained by numerical integratian. 

By dividing, point by point, the assumed deflection curve by 

this value and multiplying by the known value κ. a value of m 

for each point is obtained. If all values of m agree, the correct 


deflection curve was assumed. If they do not agree, the curve 
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yk*/m is assumed to be the deflection curve, and the process is 
repeated until all values of m are the same. Since each assumed 
deflection curve is knoun to be more accurate than the preceding 
one (6), the Vianello process uill eventually converge to a 
stable value of m. 

The Rayleigh method makes direct use of the fact that a 
column is in a state of neutral equilibrium uhen it is subjected 
to the buckling load. This can be expressed by the condition 
T D T is the work done by external loads as the column is 
displaced from a straight to a bent configuration and V ís the 
strain energy of bending induced in the column as it is displaced. 
Additionally, T = Pv, where v is the decrease in distance between 
the ends of the column as the column assumes the deflection curve 


ye For small a v is approximated by: 


coy dx, and so 
da Pf SX 


V may be expressed ЖЫ? dx where M is the bending moment. 
Lo 


ipee M = Ру, 


#£ 
on న ន / ox 
V = 4 ន ES 


Equating T and V ang slm ασ, 


pit Ax 
- У оҳ 
ЄЛ, 


4 
Ы N 
77 PR ie иб (2.8) 
x 


1% 


The solution for m consists of assuming a deflection curve y(x) 
which satisfies the boundary conditions, equations (2.1) and 
(2.2), and integrating to find N and D. It is generally 
impractical to evaluate D analytically; a numerical integration 
procedure is more convenient. An accurate value of m can be 
obtained by applying this procedure to a number of deflection 
curves. Uf these, the correct or most accurate deflection curve 
yields the minimum value of m. 

The main advantages of this system are speed and simplicity 
of Calculation. The choice of a family of deflection curves 
which can be varied to produce the correct deflection curve is 
difficult, however, and becomes even more so with the intro- 
duction of the non-symmetric column. 

Ritz provided an extension of the Rayleigh method which 
eliminates the previous objection. He employed equation (2.8) 
in devising a method to adjust systematically the parameters in 
an assumed deflection curve to yield a minimum value of т. 

The method, however, employs a comparatively large number of 
separate numerical integrations which would become quite cumber- 
some in solutions for m in non-symmetric columns. 

The investigation of the four methods was not intended to 
be exhaustive. Since each method will.yield a solution, the 
selection of the method to be used depended only upon the 
relative adaptability of the methods to the present study. The 
Vianello method with its automatic convergence was judged to be 
most suitable to the present study and is particularly adaptable 


to non-symmetric columns. The procedure outlined by Newmark 
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and the refinements which he introduced for improved accuracy were 


used in the solutions. These refinements are presented in 


Appendix A. 





Ja ២ ion Гог the Buckling Factor 

The family of columns studied is an extension of that con- 
sidered by Dinnik. The form is shown in Fig. 1. The moment 
of inertia of crgss-sectional area at distance x from the left- 


hand end is given by: 


z n 
2 = 22-0 44 О = Х/ = 9/4 


р "^0 9/4 9 xt এ (9 * alt) 


^ 
А (05е z 


"^ Ll FE (9/4 + aft) < x/“ < / O 


Dinnik restricted his studv to columns πο about midlength, 
i.e., b=e, 9=d, 2 =£/2. In the present ву, non-symmetric 
cOlumns are also included. The shape of an individual column is 
Uniquely determined by specifying 1,/Lo, n, alí , and CA. It 


is easily shown that the other parameters are given by: 


с/2 - 2/20 


| 


9/6 


d, 


2222 7? 


¿Z = (777) 9/2 


6 6 (FR) d/ f 
uhere K «(Z7 )i^ 


It is of interest to note that the postulated variation of 
I is attainable through linear taper for n=1, 2, 3, and +. 


Thus: п = l corresponds to a column of rectangular cross-section 
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Арпіс ләрип ишптол [ ০097৭ 


$4 


| 
p—— P1 κκ 
1 29 c3 


p 


with constant depth and linearly varying width; n = 2 is approxi- 
mated by a built-up column consisting of four angles connected 
by diagonals; n = 4 corresponds to a rectangular cross-section 
oolumn with constant width and linearly varying depth; andn = 4 
is represented by a solid truncated cone or pyramid. 

since the Vianello method automatically converges, the 
assumed shape of the initial deflection curve is not critical. 
However, the choice deserves some attention because the number 
of iterations required depends upon the accuracy of the starting 
assumption. The boundary conditions for a hinge-ended column 
require zero deflection at the two ends. When 1,/1, is unity, 
the tapered column degenerates into a prismatic column with a 
Sinusoidal deflection curve. From preliminary calculations and 
observations of the symmetric columns, it was noted that as 
11/1, becomes quite small, the deflection curve becomes approxi- 
mately parabolic. With these factors in mind, the initial 


deflection curve is assumed to be: 
y" D E T O A) Gl 


where u = 1 - 2x /£. Since the amount of variation in the 
deflection curve caused by the introduction of a non-symmetric 
column was unknown, the variation is ignored in equation (3.1). 
The general solution used to find the critical buckling 
factor can be best explained by an example, Table I. A typical 
non-symmetric column is chosen with n = 2, c/ = 0.25, a/ = 0.0, 
and 11/1, = 0.5. For this example, the length is divided into 


8 segments, i.e., А = 0.125 and k = 8. 
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The computed values of 1/1 are entered on Line 2. Lines 
3, 4, and 5 are self-explanatory. Line 6 is the assumed value 
of the deflection curve of the column when subjected to the 


buckling load. Line 7 was obtained by employing equation (2.6). 


y; a ($) у: (2.6) 


Line 7 is the negative of the quotient of the corresponding 
values on Line 6 by those on Line 2. Line 8 is the result of 
altering Line 7 with a refinement by Newmark. This refinement 
is explained in Appendix A. In a symmetric column, the slope 
(y!) is zero at midlength. For a non-symmetric column, this 
condition is not met, but it is nevertheless used to obtain an 
initial value for Line 9. Using this assumption, the initial 
value on Line 9 at x = 0.1582 is set equal to one half the 
negative of the entry at x = 0.5.0 on Line 8. The entries to 
the left of x = 0.55.6 are determined from right to left by 
successive subtraction of the intervening numbers of Line 8. 


Line 10 is determined by employing equation (2.5). 


уг = у,” λγε-. (2. 


The initial value is obtained from the boundary condition that 
y =O at x = Q. The following entries are the cumulative sum 
of the entries of Line 9. The incorrect choice of the position 
of zero slope and the incorrect assumption of a symmetric 
deflection curve lead to a non zero entry at x = £ on Line 10. 


The second boundary condition specifies that this value be zero. 
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Line 11 provides a linear correction, tal 106 ΕΠ“ ΠΗ. ие 
boundary condition to be met. The entries of Line 11 are com- 
puted by multiplying the negative of the x = £ entry of Line 
10 by the respective entries of Line 1. The entries of Line 
12 are the sum of the corresponding values of Lines 10 and 11. 
Finally, the values of m on Line 13 are the quotient of the 
corresponding entries on Line 6 by those on Line 12 multiplied 
by к“ (ЕСЕ 

since it is desired to determine m to four significant 
figures (three decimal places), the criterion for convergence 
is that the maximum and minimum values of m on Line 13 be 
Within 0.0005 of each other. If, as in Table I, this has not 
Occurred, Line 12 becomes the neu assumed deflection curve 
(Line 6) and the process is Tepeated to convergence. At con- 
vergence, the value of m is taken to be the average of the 
maximum and minimum entries of Line 13. 

This general method of solution was programmed for the 
COC 1604 high speed digital computer. The factors a/f , 
ΞΕ. ТИТ. and n were varied to cover the usable range 
Br COlumms of this type, amd thegresulting ous Ling factors 
are listed in Table II. The factor n varies from 1 to b in 
౭౫ steps. For each n, 6/2 varies from 0.00 to 0.50 in 
increments of 0.05; it should be noted that values of c/Z£ 
From 0.5 to 1.0 are merely repetitions of those from 0.5 to 
ЕО. гог cach r P MM 1,/1, Varies fron deel ౨౦౦ in incre- 
ments of 0.05 and a/# varies from 0.0 to 2c/@ in increments 
OF O.1. It επαω ΠΡ Егера аа ит យ defined when 


a/2 4 is greater than c/. 
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b. Discussion of Results 

ln a preliminary investigation the factor k was varied 
maca in the program solution for a column which 
appears to have one of the most stringent requirements for | 
correct convergence. It was assumed that an adequate value 
ОГ К for this column is adequate for the remainder of the 
CULL hescolunmeis Symmetric, 1.8% c/f = 0.5; and 


п = 11/] = 0.10, and a/f = D.D. The following converged 


2 


values af m were obtained for each k: 


a) k = 20 m= 4.81242 
ED) NN п = 4.81206 
πι m = 4.81203 
d) k = 160 т = 4.61203 


It is evident for the six significant figures shown, that when 

k = 80, the process converges to the proper value. The require- 
ment is for only four significant figures, however, and k = 40 
Safely fulfills this requirement. 

The final results in Table II were analyzed by: verifying 
that the solutions for prismatic columns, i.e., a// = 1.0 ог 
11/1, = 1.0, agree with theory; verifying that the results agree 
With the few analytical solutions for which there is a known 
deflection Curve; and comparing the results with those buckling 
Factors available from Dinnik's solutions. 

For a uniform slender column with hinge supported ends, the 
buckling factor m is equal to Ж. The buckling factors in 
Table II agree with the value of т“ to four sionificant figures 


For all columns where either 2/4 ӨТ 11/1, Ви = 
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There are four analytical solutions which are both per- 
tinent to this study and are readily obtainable: 
C1) αἱ =0.0, πε, NL UNUM D=1,/I,<1.0 


౧/౮ 70553 қ 


2, af -0.0, 0 ж1,/1,<1.0 


τ. σρ -П.П, пен; ЕН 0 *1,/1,+<1.0 
(КУ αἱ =0.5, πε, তা Пе 1,/1,<1.П 

These solutions were adapted from similar solutions by 
Bleich (12) and Timoshenko and Gere (13) and are presented in 
Appendix В. A comparison of the buckling factors obtained 
From them with those of identical columns from the numerical 
solution in this study iS contained in Table III. The agree- 
ment between the two is excellent. The two non-symmetric 
Cases have a maximum deviation between analytical and numerical 
results af 0.00002. This is quite good since the criterion for 
convergence in the numerical program was that the maximum and 
minimum values be within 0.0005 of each other. The maximum 
deviation of П.ПООІП in the two symmetric cases is not as small, 
but is still well within an acceptable limit. 

The buckling factors which are applicable from Dinnik's 
work were obtained from Timoshenko and Gere (13) and are 
listed in Table IV. Under each of these values is the corres- 
ponding value from this study. [here is a small deviation in 
these entries in a large number of cases; n - l seems to be the 
prime offender. The maximum difference is 0.06. Since the 
present study agrees so closely with the exact solutions, and 
Since the same method of solution is employed at all times, 
Dinnik's results seem to be in error by the amount of deviations 


from the present study. 
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ANALYTICAL CHECK UF THE CRITICAL BUCKLING FACTOR 
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2/4 


COMPARISON WITH DINNIK'S VALUES 
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‘oper values in each set of numbers are from Dinnik's work; 
those from the present work are directly below. 
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The results, as mentioned previously, cover the usable 
range for columns of the type under study. There might be a 
necessity, however, for obtaining intermediate values within 
this range Which are not listed in Table II. Therefore, some 
consideration should be given to the method of interpolation. 
Since the data obtained form curves of second order or higher, 
an error is introduced when employing linear interpolation. 
The amount of error depends upon the position of the needed 
value in Table II. The error increases with increasing n, 
decreasing a/f, and decreasing 11/12* 

A comparison of linear interpolation with interpolation 
using Bessel's formula is given in Appendix D. Interpolation 
formulas and functions were obtained from Comrie (14). For a 
μμ C =045រ ក = យ 1,/15 = 2.25, апа a% 0.15. 00 
uas found that there is a 0.11% difference between linear and 


second degree interpolation. 
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APPENDIX A 
NEWMARK REFINEMENTS 
Newmark (6) presents refinements which improve the accuracy 

of the Vianello method. He bases these refinements on the 
assumption that y" is a second degree curve, whereas, in the 
basic method, it is assumed to be linear. He assumes that 

this curve takes the form of a distributed loading curve, Fig. 
2(a), and then solves for the reactions at the end of each 


inerement (A). He finds that: 


A, = 2 (79 + бг -5) (A.1) 
πο. 24 (33 + /OF -35) 
πο τ A 3s ον 
ο ο ο. 


A. iR (Ts *Ór - 9) 


Equation (A.2) is used in the present study uhen the curve 
is smooth. When there is a break in the curve, equation (A.1) 
is used. The breaks in the curve appear where the tapered 
sections meet the prismatic section or each other. The curve 
is assumed to take the shape of Fig. 2(b) and the modification 
18: 
Ri = (749 morgens - ৩) 


The reaction R is the concentrated angle change (A y!). 


ps 








Figure 2 Newmark Reactions 





APPENDIX B 
ANALYTICAL SOLUTIONS 
Че =ОШ ЕТ Е ЕШ, Пп - 2, а/) - П.П, 0.041,/1,<1.0, 2/2 BE 
The column is of the form shown in Fig. 3, and: 
Leh (F) --ᾱ δ᾽ 
Ж -2\ / Е А Š 
Ihe differential equation of the deflection curve is: 


σ/ 


веет вета. For 1. 





& g^y 2 
2 

where ex = — 0 

22 
The general solution of equation (8.1) with two arbitrary constants 
A and B is: 

P 7 LA s/n (5 ο Θες "រ! (B.2) 

where 


А = 27 


The boundary conditions for a symmetric hinged end column may be 


PA 
;- (z 
= / ) %- O at AZ 4/2 


Substituting these boundary conditions into equation (8.2) 


taken as: 





y LO оа х= О 


results in two homogeneous equations for the constants A and B: 
+ 29% Б 
A sinlh /7/27)5 1৮8 ০০০1৮ 471%/+ |=០ 
A Z 


Ah + B⁄2 = O 





80 





0*0=4/8 
001০2 ০1790001465 ç 91π3τ1 


N 





SI 


solutions for A and B exist only if: 
tan|h (Е) | -2h=0 
2 


| ౫౯౩11 ౯౧ mumber of roots h. The smallest root ul 





defines the critical buckling load. This can be written as: 


А 4 


yo ΟΛ 


where 
πα, 


With a given ratio of I/I», Y can be obtained and the root n is: 


р = 2৮ 
фа) 


апа x NUS 
ОО 
(b) Non-Symmetric Column, n = 2, aA =0.0, D.0s [1/12- ES 
យ От 1.0. 
The column is as shown in Fig. 4 and, as for the preceding 


column, the general solution is: 
7 =, [Asin lh In$) + Bos (bln €)| 
The boundary conditions are: 
DAÑE 
5 -(%; y = Q at Xx 


រ: у= 4t X 
These conditions yield the two equations: 


A sinlh In (1/,)*] +B cos [h Ln (LM,)*]=0 


ΞΟ 
Therefore, 


sin[h {η (1/1) ] о 


|| 
O 


॥ 
N 
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С 


0°1=4/° “0*0=47е 
ишптор S511i4əmmKS-uoN ἢ 31N3T Y 
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ОГ 


ЕЕС. 
h, = In C4, 7L) 


же = |! - (1,/1,)* ΤΕ Lia 


(c) ១5៣ 7 итп, п = 5, а/Й = П.П, c/ =0.5, 0.0= ү/1„* త 


The column is of the form known in Fig. 3 and: 
*\ 4 
Е x 
Ir = 7,(5 


The differential equation is: 


El, The + Py = 


substituting = VAR 


22 
| ТҰ я | и ДК 


2 
Y 
ЕЕ +e =o (8.3) 
where „4 
_ pi? 
$^ ក 


Equation (8.3) is a form of Bessel's differential equation and 


has the solution: 


2 


s 


[4,080 + B, Y- ఊదు] 


L 
2 
where: 


Cos (SO 
J-4 (80 T TAREA 


р _ రధం 
1 ౮౦ = శరణ 
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y = ESFs |A cos (60) + Bor (St) 


and the general solution is: 


y = x” A. cos(S) + Bosin ) πο) 
The boyndary conditions are: 
“=, y = Оо Ж x 2/6 
- 5 ο’ 
a  , s= o at X = L2 


Substituting these boundary conditions into Equation (8.4) and 


selving for Ар and B yields: 


telt) = F 


or 
tan B এ 
/8 í 
where 
= = #12 
Then, 


- 448% 2,2% 


(0) Non-Symmetric Column, n = 4, a/f =0.0, c/@ = П.П аг 1.0, 


< = 
О.О= 1/1, រិ 
The column is of the farm shaun.in Fig. ^ and the general 


solution is the same as for the preceding column 


y = "A, cos(52) + В, 2%) সর 


65 


The boundary conditions are: 


χ- -, y 


8 27 х= ο 


o 


"UR T y 


Ó 27  »% =. 


Substituting these boundary conditions into Equation (B.4) and 


Solving for Ao and B5 yields: 


(84) - 


which has an infinite number of roots 77. The smallest non- 


trivial root 7 defines the critical buckling load. Then, 


D Вз 
By 


and 


А 7° 
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APPENDIX D 
INTERPOLATION METHODS 
A value of m is desired for c/( - U.5, mn sd m In R 
ఆ 
(a) Second degree interpolation 


Bessel's formula is modified to: 


m = m pAg + BEL ( вл”) 


О 


where it is assumed that the third differences are negligible. 


ЕЕЕ а 


ళీ $ 


a/l m A A 
គ 5.491 501) 
ក. 67085 654 60 


Nn EE 
D.3 7.452 
95. 085. = 0.5 СББ Ша ав 


693 2 


I,/I, = 0.20, af = 0.15 


// 
a/f i ΛΙ 
౦6 నా. 
ប្រ! 6.589 608 39 


0.2 Ja 16 
Wa 7.821 


π 589 . [5 (.608) = αρ ΘΙ 


624 


11/1. - 1,25: ЕШ ΠΡ 


a/f m N N 
ព 6.461 536 

E 6:39 7 560 24 
ШІ? y 562 2 


a2 Oe ule 
m = 6.997 + 0,5 560) = ШЕИ ШАВИ ЕУ ЕСУ 


D 


I,/I, - 0.30, a/f = 0.15 


/ ai 


8/4 m A A 
ШЕШ x 500 

Oak on 511 14 
Вр 7.095 505 -9 


0.3 επι 
7 ) 7! 0.5 (Celi) 2.062 (085) = 7.5398 


1,/1, = 0.23, a/f = 0.15 


1 1 /1 2 m А. A” 
0.15 6.406 LBL 

0.20 6.8990 385 -99 
៣. 2275 393 -62 


0.30 7.598 
m 6.825៣ (2215) 2 050 т 


(b) Linear Interpolation 
7 তা 
m = 6.589 + „304 = 6.893 


11/1, 


m = 6.997 + Olm (.560) - 7? 


= 8.28, 8/4 = BS qus 


1,/1, = 0.23, &/f = 0.15 


m = 6.893 + 0.6 (.384) = 7.123 


(c) Comparison between the two methods 


Е _ πο - 2.1642 — 
difference - — = 911% 
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